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Meth od of. and Software for. Conducting Motion Corr ection for a 

T'/tiMA^apliip Spanner 



Field of the Invention 

This invention relates to a method of, and computer software for, 
conducting motion correction for a tomographic scanner, in particular but not 
exclusively a positron emission tomography (PET) scanner. 

BackgTQiiTid of the Invention 

A typical emission scan using a PET scanner begins with the injection of 
a solution including a tracer, which is a pharmaceutical compound including a 
radio-isotope with a short half-life, into the subject. The subject may be human 
or animal. The tracer moves to, and is typically taken up, in one or more organs 
in the subject according to biological and biochemical processes which occur 
within the subject. When the radio-isotope decays, it emits a positron, which 
travels a short distance before annihilating with an electron. This annihilation 
produces two hi^ energy photons propagating in opposite directions. The PET 
scanner includes a photon detector array airanged (usually in a ring-shaped 
pattern) around the scanning area. If two photons are detected within a short 
timing window, a so-called "coincidence" is recorded along a line of response 
(LOR) connecting the two detectors. Coincidence counts , along each LOR are 
listed each time a coincidence is detected between the corresponding detector 
pair, and subsequently the Ust is processed by, for each event, incrementing a 
data storage part referred to as a sinogram bin to build up a sinograna. The 
output sinogram is typically processed using image reconstruction algorithms to 
obtain volxmietric medical images of the subject. 

For quantitative results from PET images, attenuation correction 
generally forms one of the data correction stages. In a conventional PET 
scanner, the scanner is provided with one or more positron emitter rod sources, 
formed of a material such as **Ge, which emit dual annihilation photons. To 
derive attenuation factors, two acquisitions using the rod sources are 
conventionally used - a blank scan, in which the subject being scanned is not 



present in flie scanning area (typically, the scanner is empty excqpt for the 
presence of the sources) and a transmission scan in which the subject is present 
in the scanning area. Since the source material is a positron emitter, the two 
photons arising from the annihilation of a positron and an electron are acquired 
5 in coincidence, in tihe same manner as with an emission scan. Conventionally, 
the results of the blank scan are then divided by the results of the transmission 
scan to derive an attenuation sinogram. The attoiuation sinogram is then used to 
correct tiie emission scan for attenuation. 

Typical PET scanners have detector arrangements which do not rotate 

10 during an acquisition. Typically, the detectors are arranged in,two or more banks 
of detector rings. Alternatively, the detectors may be arranged in a non-ring- 
shaped pattern. In any case, there will generally be directions in which 
coincidences are not detectable due to the geometry of the detector array, since 
the scanner has a finite field of view and there may be 'l^lind spots" due to gaps 

15 between the detectors, etc. 

With the increasing resolution of PET scanners, subject movement 
becomes an important degrading factor in the quahty of the data. For example 
with brain scans, head movement causes a time-varying rigid body 
transformation of the brain, and hence of the radioactivity distribution to be 
reconstructed. The head movement can be extemally monitored, and the effects 
of the monitored movements can be corrected Several methods have been 
proposed to correct for subject movement. They fall into 3 classes: 

(a) Post-processing of the image using deconvolution. See for 
example the papery by M. Menke, M.S. Atkins, K.IL Buckley, "Compensation 
Methods for Head Motion Detected During PET Imaging", IEEE TNS vol. 43, 
(1996)310. 

(b) Multiple Acquisition Frame (MAF) methods, i.e. spUtting the 
acquisition into short time^fi-ames (which may be triggered by a threshold on the 
amount of motion), followed by reconstruction each short time fi-ame, 
realignment of the images and then combining all these images into one image. 
See for example the paper by Y. Pickard, C.J. Thompson, ^Motion correction of 

« 

-PETHbtnages nsing^mxdtiple acquisition fram^^ IEEE TMrvol.~l'6(1997)"1377 



(c) Acquiring the data in list-mode and realigning the events 
according to the known motion before binning them into a sinogram, followed 
by reconstruction of the sinogram. See for instance the paper by P. Bloomfield, 
TJ. Spioks, J. Reed et al., "the design an4 implementation of a motion 
correction scheme for neurological PET", PMB 48 (2003) 959 for a discussion of 
the merits of these methods and more complete references, and a paper by S.K. 
Woo, H. Watabe, et al, "Sinogram-based motion correction of PET Images using 
Optical Motion Tracking and List-mode Data", Conf. Proc. IEEE MIC 2002. 
Both the Bloomfield et al. paper and the Woo et al. paper describe reconstructing 
the data with a Fourier Rebinning (FORE) algorithm followed by a Filtered Back 
Projection (FBP) algoriliim, and mention that artefacts in the images depending 

on the type of motion. 

A paper by H. Watabe, N. Sato, et al, "Correction of Head Movement 
Using Optical Motion Tracking System during PET Study with Rhesus 
Monkey", "Brain Imaging Using PET" eds. M. Senda, Y. Kimura, P. 
Herscovitch, Academic Press ISBN 0-12-636651-9 pl-8 2002 observes tiiat in 
2D PET, rotation can move LORs to oblique LORs, and suggested a solution 
which essentially amounts to a Single SMce Re-Binning (SSRB) approximation. 

It is an object of the present invention to provide improvements relating 
to motion correction for tomography, in particular but not exclusively PET 
scaxmers. 

In accordance with one aspect of the present invention there is provided 
method of conducting motion correction for a tomogr^hic scanner including a 
detector array for detecting radiation to generate detector data, wherein the 

method comprises: 

storing detector data collected .during a data acquisition period, said 

detector data being iiidicative of: 

i) directions along which radiation is' detected; and 

ii) quantities ofTadiation detected in different of said directions; 
Storing movement data representing movement of the subject durmg the 

data acquisition period; and 
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motion correcting said detector data using said movement data and a 
motion correction algorithm to calculate motion corrected detector data, 

wherein said motion correcting step comprises processing said detector 
data by: 

a) realigning directions of at least some of said detector data on 
the basis of said movement data; and 

b) altering quantities of at least some of said detector data on the 
basis of said movement data, 

such that at least some of said detector data are both realigned and altered 
in quantity. 

In preferred embodiments of the invention, the altering quantities step 
comprises scaling quantities upwards and/or replacing quantities with calculated 
quantities. Prefffably, the method comprises selectively either scaling a quantity 
upwards or replacing a quantity with a calculated quantity. Altematively, the 
method may comprise altering the quantity so that it the altered quantity takes 
into account both the original quantity and a separately calculated quantity (for 
example, by weighted average of a rescaled quantity and a separately calculated 
quantity.) 

Typically, the motion corrected data will then be used as input data for an 
image reconstruction algorithm to produce a reconstructed image. The present 
mvention provides a new motion correction procedure which can be used to 
reduce the appearance of artefacts in the reconstructed images due to particular 
types of motion of the subject, including translation along the scaimer axis and 
rotation about a rotation axis orthogonal to the scanner axis. 

Aspects of the invention further include computer software and a scanner 
system arranged to carry out the method of the invention. 

Further features and advantages of the present invention will become 
apparent from the following description of preferred embodiments of the present 
invration, made by way of example only with reference to the accompanying 
drawings. 



1 



Brief Description of Drawings 

Figure 1 is a schematic diagram of a PET scanning facility; 

Figure 2 is a schematic diagram of data processing components in a PET 
scanner and an associated data processing terminal; 
5 Figure 3 is a schematic cross section of the arrangement within a PET 

scanner during an emission scan acquisition; 

Figure 4 is a flow diagram illustrating a general data processing 
procedure conducted, including motion correction according to an embodiment 
of the invention; 

10 Figure 5 is an illustration of areas of missing data in a sinogram due to 

axial linear motion of the subject during an ^nission scan acquisition; and 

Figure 6 is an illustration of areas of missing data in a sinogram due to 
non-axial rotational motion of the subject during an emission scan acquisition. 
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15 Detailed Description of Preferred Embodiments 

Referring now to Figure 1, a PET scanning facihty, arranged in 

accordance with an embodiment of the invention, is shown, which includes a ,3* 

PET scaimer 2, a tracer generator module 4 and an operator computer terminal 6. 

The scanner 2 includes a detector array 8 arranged about a scanning area, in 
20 which a subject 10 is located during a transmission scan and during an emission 

scan. 

One embodiment of the invention relates to a non-rotating PET scaimer. 

For example, a PET scanner of the ECAT EXACT3D™ type is used. In tiiis 

example, the detectors in Ihe detector array 8 are arranged in square detector 
25 blocks, each containing multiple detector elements. The detectors are arranged 

in multiple rings of detector blocks, the rings being arranged adjacent one 

another along a scanner axis. 

As shown in Figure 2, the PET scanner 2 includes a contrbl mut 30, 

detection djata processing circuitry 32, a motion detector 34, a data store 36 for 
30 storing the detector data and movement data, and an Input/Output (1/6) device 

38. 



In this embodim^t, the motion detector 34 is in the form of a Polaris 
measurement device (produced by Northern Digital Inc. of Waterloo, Canada), 
which measures the motion of a plate with foxir reflectors using infrared 
radiation. This plate is attached to a neoprene cap closely fitting the head of the 
subject, so as to track movemCTit of the head, to produce movement parameters 
indicating translational motion along three ojlhogonal axes and rotational motion 
about three orthogonal axes. 

The computer temiinal 6 includes a central processing imit (OPU) 42, 
memory 44, hard disc drive 46 and I/O device 40, which facilitates 
interconnection of the computer 6 with the PET scanner 2. Operating system 
programs 48 are stored on the hard disc drive 46, and control, in a known 
maimer, low level operation of the computer terminal 6. Program files and data 
50 are also stored on the hard disc drive 46, and control, in a known manner, 
outputs to an operator via associated devices. The associated devices include a 
display, a pointing device and keyboard (not shown), which receive input from, 
and output information to, the operator via fiuther I/O devices (not shown). 
Included in the program files SO stored on the hard drive 46 are a motion 
correction software application 52 and an image reconstruction software 
application 54. A database 56 is used to store the detector data and movement 
data transferred from the PET scanner 2. 

In acquisition mode, the detection data processing circmtry 32 processes 
all events dejected in the detector array 8, and by using a coincidence timing 
window, detects coincidences between the events which are recorded as real 
coincidence counts. These coincidence coimts are then output in Ust-mode to 
data store 36, where it is stored a Ust-mode file for subsequmt processing. In 
parallel with the data acquisition by the radiation detectors, position data is 
collected at regular, sub-second intervals (e.g. lOHz) from motion detector 34 

>and also sent to data store 36 for subsequent processing. ^ 

Figure 3 illustrates features of operation of the PET scann^ during an 
emission scan. During an emission scan, the subject 10 is placed in the scanning 
area and contains the tracer gen^ated by the tracer generator module 4. Figure 3 

"~shows"a"positroiremissionnevent being-registered in- therdetector ■array-S.—Tlie— • 



positron 60 annihilates and generates a first photon 62 travelling in one direction 
and a second photon 64 travelling in an opposite direction. The first photon is 
detected by a detector element 70 on one side of the detector array 8, and the 
second photon 64 is detected in a different detector element 72 on the other side 
of the detector array 8. As described above, the detection data processing 
circuitry 32 registers the two as a coincidence along the LOR defined between 
the two diflfrarent detectors 70, 72. 

The motion correction procedure carried out using motion correction 
software apphcation 52 includes, in one embodiment of the invention, the stages 
illustrated in Figure 4. These are the stages of acquiring the detector data in list- 
mode 100, weighting each list-mode event with the normalisation factor 
appropriate for the original detector pair 102, performing motion correction 104, 
binning the realigned list-mode events into sinograms 106, and image 

reconstruction 108. 

In the realignment during motion correction 104, for each event in the 
list-mode file, the LOR is repositioned according to the motion parameters 
derived from the motion detector 34. This can be done at sub-second time- 
resolution, since the list-mode data has a time resolution which is higher than 
that of the motion detector data (which is itself available at frequencies up to 
20Hz.) The reoriented LOR is then stored in an ^propriate sinogram bin using 
an interpolation method such as nearest neighbour interpolation. 

Since the scanner has a finite field of view the motion corrected LORs 
can sometimes fail to fill the whole 3D sinogram, with parts of the sinogram 
being completely missing or having too low values because they are 'filled' only 

partly during the whole acquisition. 

Two particularly problematic cases have been identified: axial translation 
and rotation around an axis which is not parallel to the scanner axis (for instance 
caused by the subject nodding). These two problems are illustrated in Figmes 5 
and 6 respectively, which show schematicaHy a uniform cyhnder phantom 200 
placed in the scanner, and moved between two discrete positions, from a first 
position 200A to a second position 200B (Figure 5) and a third position 200C 
(Figure 6). 



Translation along the axis of the scanner 

Figure 5 illustrates the case in which the subject has moved outward (i.e. 
.away froni.the bed end. furthest firom.the scanner), and the motion correction is 
such as to tend to reposition it towards the centre of the scanner. There is no 
measured data that will be repositioned into flie set of non-obUque LORs at tiie 
end of the detector array. So, the reconstructed image will have problems in the 
end-planes. 

Further, in a 3D PET acquisition,, data for oblique LORs are also 
acquired. Obhque LORs are those in which each end of the LOR is on a 
different ring (the number of rings separating the respective sides of the LOR is 
referred to as a ring-difiference.) For each ring-diflference, the same number of 
obhque end-LORs will be affected, as can be seen on Fig. 5, which illustrates the 
missing oblique data for a selected ring difference. These missing obhque end- 
LORs would have, had the subject not moved, measured parts throughout the 
subject which are not near the ends of the scaimer. In fact, for the maYi'mnm ring 
difference, the missing LORs go right through the centre of the scanner axis. If 
the fact that these end-LORs are erroneous is not taken into account, the 
resulting image will show artefacts. The type of artefact dqpends on the 
reconstruction algorithm used, as each algorithm reacts differently to 
inconsistent data. 

...As an example, a 'sagital' shoe was taken through the reconstructed 
image of a motion corrected cylinder, i.e. a shoe through the axis of the scanner. 
An analytic reconstruction algorithm called three dimensional reprojection 
(3DRP) was used. Artefacts were clearly visible in the object relatively far from 
flie end-planes. This can be explained by the edge enhancement of the Colsher 
filter togethCT with the fact that, in this expmment, the data for the end-LORs for 
a fixed ring difference were underestimated, while the repixgected data were not 
The artefacts obtained using the FORE analytic reconstruction algorithm were 
less severe, but were still clearly visible in areas away from the edge of the 
cylinder. 
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Rotation around an axis orthogonal to the scanner axis 

Figure 6 illustrates the case with rotation of the subject about an axis 
orthogonal to the scanner axis. On reaUgnment during the motion correction 

5 procedure, the set of oblique LORs, those making angles relative to the line 
connectrag diametrically opposite detectors in the plane of the moment within 
the angular extent of the rotational movement, will have missing data. 
Therefore, if the rotation is too large (i.e. larger than the axial acceptance angle 
of the scanner), the motion corrected segment will again cause problems in the 

10 throughout the scanner, not only in the end-planes. 

Continuous motion 

In this case, the sinogram resulting after the binning process can be seen 
to be simply a superposition of sinograms obtained for very small time duration. 
15 Some bins in the projection data produced by the rebinning procedure described 
above will have the same value as if no motion had occurred, but others will 
have a smaller valxie, depending on how many of these very short time frames 

would fill data in the bin. 

In the notation used herein, elements of the (3D) sinogram are indexed by 

20 b (for bin). The projection data that is obtained for the subject in the position at 
time t will be denoted by s\ and the reference position, at which the reference 
sinogram is taken, is denoted with a subscript r. The number of very short time 
frames is denoted by The motion correction procedure of binning each Ust- 
mode event into a sinogram bin can be seen as a linear process which can be 

25 conceptually encoded in a matrix that will change with t: 

s^I^Ym's' (1) 

Note that this matrix M' combines the repositioning of the LORs with the 
binning into a sinogram. So, it would have dimensions nxw where n is tiie 
number of sinogram elements after binning, and m is the number of possible 
30 LORs, with n<m being the usual case. 
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Each motion corrected sinogram shoxild be the same as the reference 
sinogram, except for Ihe missing data, at least if discretisation and interpolation 
errors are ignored. So; 

M^s' ^diag(d^)s' (2) 

where diag(d^) is a diagonal matrix with elements d'b on the diagonal which are 
0 where the data are missmg, 1 for bins which were completely in the field of 
view, and potentially between 0 and 1 for bins at the border between these two 
regions. So, the motion corrected sinogram will have elements which are scaled 
versions of the reference sinogram: 

Sb=^d^5''b9 (3) 

with: 

^"^j^lL^'o (4) 

t 

Since these scale factors are not all equal to each other, the motion 
corrected sinogram will be inconsistent (i.e. does not correspond to the sinogram 
of the subject in its reference position) and the reconstruction process will have 
to be adapted to take this into accoimt. Otherwise, artefacts will occur. 

Making the motion corrected sinogram consistent 

When a scale factor is non-zero, we could just divide the detector data 
sj, by it. However, for small scale factors division would lead to noise 
enhancement in the data. Therefore, a threshold is used so that when the scale 
factor is smaller than a threshold, the corresponding bin cannot be corrected by 
rescaling. Instead, an estimate of the value which should be in the bin is 
separately calcxilated. From an initial estimate of the image, taken from only the 
non-obUque LORs, referred to as the direct sinograms. This is then forward 

projected to obtain an estimate of the value s^b of the bia in the reference 
sinogram. 

A consistent sinogram can be estimated as: 



2 J i:aii.1'X 

where s-^ t is the value of the forward projection of the initial image, and uj, is 

the threshold The threshold can be a global threshold appUed to all bins, or may 
be a local threshold varying depending on the bin identity. 

The motion correction software apphcation 52 preferably includes 
5 functionality allowing the threshold, whether global or local, to be varied 
between different acquisitions, depending on the quality of the data obtained. 

Construction of the initial estimate of the image 

As mentioned above, the initial estimated image can be produced by a 

10 reconstruction of the direct sinograms. However, even the direct sinograms can 
have missing data, especially in the end planes. One solution for this is to 
discard those raid planes of the initial image that are too corrupted. In the case 
that too many direct sinograms are corrupted by rotation, an alternative would be 
to reconstruct the initial image from a set of oblique projections orthogonal to a 

15 direction not parallel to the scanner axis. 

With both these procedures, it is possible to end up with an initial image 
that is shorter than the scanner field of view. It is then important to consider 
which LORs are used in the final reconstruction process. The values of the 
estimated bins b will correspond to an object which is as long as the initial 

20 image, while the motion corrected data Sf, could have contributions of activity 
which does not overlap with the initial image. These bins are discarded, as 
otherwise the data would stiU have potential inconsistencies. 

The end result of this procedure will be a final reconstructed image that 
will have fewer planes than usual (i.e. the same number as in the initial image), 

25 but the image will have no artefacts. 

Determination of the scale factors 

The scale factors d^h defined in equation (2) are (approximately) 
independent of the sinogram that is being motion corrected. Their role is similar 
30 to the (inverse of the) usual normalisation factors in PET. 
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The scale factors can be computed by motion correcting a data set which 
has 1 in every LOR, and a set of LORs is sampled after the motion correction to 
produce scale factors which vary across the entire sinogram. 

Both sinograms in equation (2) have now all elements equal to 1, so: 

d'=M'u^ (6) 

with a ixTO colunm vector with all elements equal to 1. This can be 

computed by a loop over all possible LORs, repositioning the LOR, and storing 
'events' in the sinogram using the same interpolation as used in the actual 
motion correction. 

Computing the scale factors J^, which are averaged over the whole 

frame using equation (4), is highly processor-intensive. Several optimisations 
are possible. The number of time steps in the sum can be made smaller as the 
motion parameters tend to vary smoothly aside from sudden jumps. In addition, 
if an appropriate interpolation scheme is used during the motion correction, the 
set of bins which has scale factor 1 is a geometrically connected region. This 
means that the scale factors can first be computed in with low precision, and ' 

that then edges in the resulting sinogram can be recomputed with higher 
precision. 

Determination of the thresholds 

The role of the thresholds is to avoid noise amphfication. A niraiber 

of approaches are possible. For instance, strong smoothing on the motion 
corrected sinogram s can be appUed to obtain a rough estimate of the noise level 
in each bin which would enable us to use lower values for the higher values in 
the each bin. 

This means they can be chosen lower where the motion corrected data is 
less noisy. Because of the Poisson nature of PET data, the relative noise level 
decreases with higher values of the mean (over repeated experiments) value of 
. However, the data in a single sinogram bin is generally too noisy to use as 

an estimate for this mean value. This can be remedied by apply a strong 
■^moothing-Qirthe-motion-correctedrsinogramr-.sH:©' obtain-a rough-estimate-of-the — 
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noise level in each bin. The threshold can then be made lower for higher values 
of the bin in the smoothed sinogram. After this, the smoothed sinogram is 
discarded. Alternatively, a global threshold (i.e. independent on b), may be used. 
The global threshold may be varied in dependence on a factor such as the total 
number of counts in a sinogram. 

The higher the threshold, the more data will be replaced by the separately 
calculated forward projections. Since the estimated forward projections 
generally have lower resolution compared to the measured data, setting the 
threshold higher will result in images with lower resolution but less noise. 



Normalisation for detection efficiencies 

In practice, different LORs will have different detection efficiencies, 
which can be used to form a detection efficiency sinogram e , The detected 
sinogram diould be normalised to correct for these efficiencies. In one 
15 embodiment this is done before binning (i.e. the detection efficiency factor is 
obtained for an event, the LOR is repositioned, the event is weighted with the 
inverse of the detection efficiency factor and interpolated into the motion correct 
sinogram). 

An alternative solution is to not perform normalisation during binning, 
20 but to include these scanner normahsation factors in the scale factors referred to 
above. In the case without motion correction, this has been shown to give better 
noise properties. To do this, the motion correction sinogram is obtained as in 
equation (1), where now is however not precorrected for detection 
efficiencies. The result will again be (almost) proportional to a reference 
25 sinogram as in equation (3). Take as reference what would have been obtained 
without motion, but in this case with normahsation for detection efficiencies. 
These scale factors can then be found by looking at smooth objects. That is, by 
assuming that the reference data is a sinogram filled with Is. The non- 
normalised measured data will then be the detection efficiency sinogram s . The 
30 analogue of equations (4) and (6) is now: 

d^^YM's (7) 



•••• 
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These scale factors will no longer all be less than or equal to 1, as they depend 
on £ . This can be taken into account when setting the thresholds. 

Other precorrections 

If the list-mode data contains delayed events, then randoms correction 
can be performed during the binning by subtracting tiie delayed events after 
repositioning. If some other way of estimating the randoms is used, the motion 
» would have to be taken into account. 

For scatter , it can be assumed that all of the scatter is due to the subject 
In diat case, the scatter distribution is just moved along with the subject resulting 
in a scatter distribution which is exactly the same as the scatter distribution of the 
subject in the reference positioru So, any existing scatter correction procedure 
can be used on the motion corrected data. However, some scatter is due to non- 
moving objects, such as the bed and the scanner itself. This can be corrected by 
estimating the scatter distribution caused by the non-moving objects in multiple 
positions and motion correcting the results. 

Finally, there is attenuation correction. However, after motion 
correction, the attenuation factors are the same as those for the subject in the 
reference position. Preferably, the reference position is chosen to be the position 
of the subject during transmission scanning. 

Description of specific procedure 

A step-by-step description of one embodiment of the motion correcting 
procedure follows: 

1. Obtain attenuation coefficients of subject in the reference 
position. This can be done using the following procedure: 

(i) "Perform a blank and transmission scan. r 

(ii) Reconstruct an image representing the attenuation map of the 
subject using the transmission and blank data. 

(iii) Move the attenuation map to the reference position by 
remteipolation. 
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(iv) Forward project this attenuation map and exponentiate. 

2. PCTform a list-mode emission scan while tracking the motion of 
the subject 

3. Bin the list-mode data , into a motion corrected sinogram by 
5 repositioning the LORs. Normalise the events for detection efificioicies. 

4. Obtain the scale factors equation (4) or (7) depending on whether 
events were normalised or not in the previous step. 

5. Obtain tiboresholds a^. 

6. Divide the motion corrected anogram by the scale factors , 

10 unless the scale factor is too small. In tiiat case, rescale the data with the 
threshold: 

si,ldf, ifdf,>af, 
Sf, / otherwise 

7. Correct the resulting sinogram s"^ for scatter and attenuation as 

usual. Call the result . 
15 8. Reconstruct initial image from a consistent subset of this 

precorrected sinogram. This image can be shorter than the field of view. 

9. Reduce tiie size of the 3D sinogram such that only LORs are 
considered which intersect the initial image, but no other parts of the subject. 

10. Merge the forward projections (by using line integrals) of the 
20 initial image into the precorrected sinogram: 

s^'b ifd^>ai, 



s'^b^ 



s'b^ 



y where ^ is the forward projected data. 
s^b otherwise 



The result of this procedure is a consistent, fully-precorrected 3D 
sinogram (ignoring the effects of noise). This can then be passed to any known 
25 3D reconstruction algorithm that works on fully precorrected data. 

Alternatives to the above-described embodiments are envisaged. Whilst 
in the above embodiment, the detector data is output to a list-mode file which is 
then processed offline to perform motion correction, in an alternative 
embodiment of the invention the motion correction system of the present 
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invention may be built into a scanner, so as to provide online motion correction 
of LORs which can then be used to directly generate a sinogram. 

Whilst in the above-described embodiments, where the original measured 
data is inconsistent due to motion of the subject, the original measured data can 
be either rescaled or replaced with a separately calculated value. Alternatively, a 
weigihted average of the two types of corrected data may be taken when the scale 
factor falls below the preset threshold; or a varying weighting can be applied in 
dependence on the scale factor value. 

It should be noted that the invention applies also to PET scanners in 
general, both non-rotating and rotating. Further, although the above-described 
embodiment relates to a PET system, the method of the invention could also be 
used in other systems of tomography, such as x-ray computer tomography (CT) 
and single photon emission computer tomography (SPECT). 

The motion correction method of the invention may be used in 
combination with analytic image reconstruction algorithms such as 3DRP, 
FAVOR, BPF, etc. Altematively, or in addition, the method of the invention 
may be used with analytic rebinning algorithm such as FORE, FORE-X, FORE- 
J, in combination with an iterative reconstruction algorithm, or indeed iterative 
algorithms which work on precorrected data 

The above embodiments are to be understood as illustrative examples of 
the invention. It is to be imderstood that any feature described in relation to any 
one embodiment may be used alone, or in combination with other features 
described, and may also be used in combination with one or more features of any 
other of the embodiments, or any combination of any other of the embodiments. 
Furthermore, equivalents and modifications not described above may also be 
employed without departing Jfrom the scope of the invention, which is defined in 
the accompanying claims. 
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1. A method of conducting motion correction for a tomographic 
scaimer including a detector array for detecting radiation to generate detector 
data, wherein the method comprises: 

storing detector data collected during a data acquisition period, said 
detector data being indicative of: 

i) directions along which radiation is detected; and 

« 

ii) quantities of radiation detected in different of said directions; 
storing movement data rq)resenting movement of the subject during the 

data acquisition period; and 

motion correcting said detector data using said movement data and a 
motion correction algorithm to calculate motion corrected detector data, 

wherein said motion correcting step comprises processing said detector 
data by: 

a) realigning directions of at least some of said detector data on 
the basis of said movement data; and 

b) altering quantities of at least some of said detector data ort the 
basis of said movement data, 

« 

such that at least some of said detector data are both realigned and altered 
in quantity. 

2. A method according to claim 1, wh^ein said altering quantities 

» 

step comprises rescaling quantities. 

3. A method according to claim 1 or 2, wherein said altering 
quantities step comprises replacing quantities with calculated quantities. - 

4. A method according to claim 1, 2 or 3, wherein the method 
comprises both rescaling a quantity and separately calculating a quantity. 
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5. A method according to claim 4, wherein the method comprises 
selectively rescaling a quantity and/or utilising a separately calculated quantity, 
on the basis of a threshold. 

5 

6. A method according to any of claims 5, wherein said threshold is 
variable. 

7. A method according to claim 5 or 6, wherein a value indicating an 
10 amount of scaling deemed appropriate for a quantity is determined, and the value 

is compared to the threshold in order to determine whether to scale the quantity 
upwards or to replace it with a calculated quantity. 

8. A method according to any preceding claim, further comprising 
15 performing image reconstraction using said improved detector data. 

9. A method according to claim 8, comprising performing image 
reconstruction using an anal3^c image reconstruction algorithm and/or an 
analytic rebinning algorithm. 

20 

10. A method according to claim 9, wherein said algorithm is one of 
.3DRP, FAVOR, BPF, FORE, FORE-X and FORE- J. 

11. A method according to any preceding claim, wherein said seamier 
25 is a positron emission tomography (PET) scanner. 

12. A method according to claim 11, wherein said detector data is 3D 
PET detector data 

30 13. Computer software for conducting motion correction for a 

tomographic scanner including a detector array for detecting radiation to 
detector 3ata, wEereinThe ¥ottwafe Is adapted'toT"* 
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store detector data collected during a data acquisition period, said 
detector data being indicative of: 

i) directions along which radiation is detected; and 

ii) quantities of radiation detected in different of said directions; 
store movement data representing movement of the subject during the 

data acquisition period; and 

motion correct said detector data using said movement data and a motion 
correction algorithm to calculate motion corrected detector data, 

wherein said motion correcting step comprises processing said detector 
data by: 

a) realigning directions of at least some of said detector data on 
the basis of said movement dat^ and 

b) altering quantities of at least some of said detector data on the 
basis of said movement data, 

such that at least some of said detector data are both realigned and altered 
in quantity. 

14. A data carrier comprising computer software according claim %3, 

15- A tomographic scanner system including a detector array for 
detecting radiation to generate detector data, wherein the scanner system is 
adapted to: 

store detector data collected during a data acquisition period, said 
detector data beiag indicative of: 

i) directions along which radiation is detected; and 

ii) quantities of radiation detected in different of said directions; 
store movement data representing movement of the subject during the 

data acquisition period; and ' 

motion correct said detector data using said movement data and a motion 
correction algorithm to calculate motion corrected detector data, 

wherein said motion correcting step comprises processing said detector 
data by: 
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a) realigning directions of at least some of said detector data on 
the basis of said movement data; and 

b) altering quantities of at least some of said detector data on the 
basis of said movement data, 

5 such that at least some of said detector data are both realigned and altered 

in quantity. 
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Abstract 

A method of conducting motion correction for a tomographic scanner 
including a detector array for detecting radiation to generate detector data, 
wherein the method comprises: ^ 

storing detector data collected during a data acquisition period, said 
detector data being indicative of: 

i) directions along which radiation is detected; and 

ii) quantities of radiation detected in different of said directions; 
storing movement data representing movement of the subject during the 

data acquisition period; and 

motion correcting said detector data using said movement data and a 
motion correction algorithm to calculate motion corrected detector data, 

wherein said motion correcting step comprises processing said detector 
data by: 

a) realigning directions of at least some of said detector data on 
the basis of said movement data; and 

b) altering quantities of at least some of said detector data on the 
basis of said mov^ent data, 

such that at least some of said detector data are both reaUgned and altered 
in quantity. 



[Fig. 1] 
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